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A suitable coupling of the inflaton <^ to a vector kinetic term gives frozen and scale invariant 
vector perturbations. We compute the cosmological perturbations ( that result from such coupling 
by taking into account the classical vector field that unavoidably gets generated at large scales 
during inflation. This generically results in a too anisotropic power spectrum of i^. Specifically, 
the anisotropy exceeds the 1% level (10% level) if inflation lasted ~ 5 e-folds (~ 50 e-folds) more 
than the minimal amount required to produce the CMB modes. This conclusion applies, among 
others, to the application of this mechanism for magnetogenesis, for anisotropic inflation, and for 
the generation of anisotropic perturbations at the end of inflation through a waterfall field coupled 
to the vector (in this case, the unavoidable contribution that we obtain is effective all throughout 
inflation, and it is independent of the waterfall field). For a tuned duration of inflation, a 1% (10%) 
anisotropy in the power spectrum corresponds to an anisotropic bispectrum which is enhanced like 
the local one in the squeezed limit, and with an effective local /nl ~ 3 (~ 30). More in general, 
a significant anisotropy of the perturbations may be a natural outcome of all models that sustain 
higher than spin fields during infiation. 



I. INTRODUCTION 

In this work we compute the power spectrum and the 
bispectrum of the cosmological curvature perturbations 
C under the assumptions that 

1. The dominant contribution to ^ is provided by a 
slowly rolling inflaton field 

2. The inflaton is coupled to the kinetic term of a vec- 
tor field so to produce a nearly scale invariant and 
frozen spectrum of vector perturbations at large 
scales. 

These assumptions are realized in models of magne- 
togenesis and of statistical anisotropy of the cosmologi- 
cal perturbations. Despite the phenomenology of these 
models has been heavily studied in the literature, we ob- 
tain novel and general results, since, for the first time to 
our knowledge, we simultaneously take into account both 
the facts that (i) a strong contribution to the anisotropy 
(both in the power spectrum and in the bispectrum of 
C) results from the same coupling characterized in 2., 
and (ii) the inflationary expansion that took place before 
the Cosmic Microwave Background (CMB) modes left 
the horizon unavoidably results in a classical background 
vector field that is homogeneous from the point of view 
of the CMB modes, but breaks isotropy. 



Interest for models that can produce vector fields dur- 
ing inflation has been generated by the inference of inter- 
galactic magnetic fields and by claims of broken statis- 
tical invariance of the CMB modes. Intergalactic mag- 
netic fields have been inferred by an apparent lack of 
GeV scale 7— rays coming from blazars that produce TeV 
scale 7— rays; in standard models, part of the higher en- 
ergy 7— rays should be converted in lower energy secon- 
daries (which then generate the lower energy 7— rays) by 
their interaction with the inter galactic medium. The 
non-observation of the GeV scale 7— rays has been ex- 
plained as intergalactic magnetic fields that defiect the 
secondaries pj. 

Broken statistical isotropy of the CMB perturbations 
has instead been found in the studies [IHl] of the WMAP 
data. While the overall WMAP results [S] strongly sup- 
port the inflationary paradigm, the above studies have 
shown that the statistics of the WMAP anisotropies does 
not possess full rotational invariance. Specifically, under 
the parametrization [6j 



P(;[k) =P{k) 1 + g, cos^ 6*^ 



(1) 



(which can be thought of as an expansion series of the 
power spectrum in the limit of small anisotropy, trun- 
cated at the quadruple term) the WMAP data give 
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= 0.29 ±0.031 [3^. The "privileged" direction V lies 
very close to the ecliptic poles. This strongly suggests a 
systematical origin of the effect, and it has been shown in 
[31 [7] that the instrument beam asymmetry can account 
for it. Fortunately, the Planck satellite will soon provide 
an independent test of this with an expected sensitivity to 
a quadrupolar anisotropy in the power spectrum as small 
as 0.5% at 1(7 [8, 9 . On different scales (and marginaliz- 
ing over the preferred direction V) Large-Scale Structure 
data analysis constrain —0.41 < 5* < 0.38 at 95% C.L. 
[lOj (the amplitude of the anisotropy may in general be 
scale dependent 

Broken rotational invariance could be the result of 
anisotropic inflation |llj . It is however nontrivial to 
realize this, since anisotropic spaces typically rapidly 
isotropize in presence of a cosmological constant [T^ . [j 
Vector fields may in principle support the anisotropy. In 
this case, the problem of preserving the anisotropy trans- 
lates into contrasting the quick decrease of the vector en- 
ergy that takes place for a minimal Ca = —F'^/A. To 
our knowledge, four distinct classes of models have been 
constructed to achieve this; the first three of them are 
characterized by (i) a vector potential ^ (A^) [13], (ii) 
a fixed vector vev due to a lagrange multiplier [5], and 
(iii) a vector coupling A^R to the scalar curvature R 
I16j . This last mechanism was originally employed 
for magnetogenesis in [T7]. These three proposals break 
the U(l) symmetry of the minimal action, and lead to 
an additional degree of freedom, the longitudinal vector 
polarization, that in all of these models turns out to be 
a ghost [IS]. 

The fourth class is instead U(l) invariant and free of 
ghost instabilities. It is characterized by a function of a 
scalar inflaton tp multiplying the vector kinetic term. 



(2) 



A suitably chosen evolution for (/) during inflation re- 
sults in a (nearly) constant vector energy density, and 
therefore in a prolonged anisotropic expansion [19]. 
Also this mechanism was originally suggested for mag- 
netogenesis psj (this application is however problematic 
[26] . as we discuss below). For anisotropic expansion, a 
homogeneous vector field pointing along a given direc- 
tion corresponds to an "electric" component, and ([2| en- 
joys an "electric" -O- "magnetic" duality under P ■>r^ -p 
[27j . A constant "electric" component is also produced 
through ([2]) in the mechanism of [55] , in which the vector 



^ See 1131 for the extension of the study of 'I2^ to slow roll inflation 
^ See |20) for models of anisotropic inflation that employ the idea of 
|19l . Also, an interesting model of vector curvaton 1211 employing 
a varying mass m and kinetic function / has been proposed in 
|22l and studied in 1231 124| . In particular, ref. 1241 demonstrated 
that treating / and m as functions of a quantum inflaton field 
results in a different phenomenology than just treating them as 
classical external functions. 



field is coupled to the waterfall field x of hybrid infla- 
tion through a ~x?A^ interaction. Due to this, the gauge 
field provides a contribution to the mass of x, concurring 
to determine the moment at which inflation ends, and 
- thanks to this - contributing to the curvature pertur- 
bation. The waterfall field acts as the medium through 
which the anisotropy in is communicated to the in- 
flation; however, as we shall see, the communication al- 
ready occurs through the very same interaction ^ that 
supports the vector field. This unavoidable effect has 
not been accounted for neither in |28| , nor in the related 
works 



The linearized theory of cosmological perturbations in 
the anisotropic inflationary model of [19 was worked 
out in [5iH57] . The classical equations of motion of the 
model admit an attractor solution, [191 El] character- 
ized by a non- vanishing "electric" component E^'^'^. A 
10% level anisotropy {\g^ \ = O (0.1)) is found for an en- 
ergy |_B*^°)p/2 which is about eight orders of magnitude 
smaller than the inflaton potential [35 37^. Therefore, 
the vector energy needs to be highly subdominant not 
too produce a too strong anisotropy. The work |39j com- 
puted instead the cosmological perturbations in the case 
in which ( 2| provides scale invariant "magnetic" compo- 
nents of the vector field, as in the magnetogenesis appli- 
cation [551 HO] ■ Cosmological applications in this context 
have also been studied in [41]. Studies of the cosmo- 
logical perturbations in [551 HO] start from the point of 
view that the statistics of the generated "magnetic" field 
is isotropic, and therefore obtain statistical isotropic re- 
sults. In shorts, in the magnetogenesis context [25l |40] 
one does not have the analogous of the attractor solution 
E'^^'^ of the classical equations of motion of the anisotropic 
inflationary model ^9j (the corresponding B'^'^^ vanishes 
in [25l |40] ) and therefore it is simply assumed that <?* = 
in this case. 

However, in [HI [28] (respectively, in [25l |40]), the 
CMB perturbations are affected by a classical "elec- 
tric" field -^classical (respectively, classical "magnetic" 
field i?ciassicai) which is in general different from the value 
given by the classical equations of motion. Indeed, such 
mechanisms are designed to result in a nearly scale in- 
variant spectrum for the "electric" (respectively, "mag- 
netic") perturbations. Let us denote by A^tot the number 
of e-folds of inflation. The modes that left the horizon in 
the first A'tot — N e-folds of inflation add up as a classi- 
cal background from the point of view of the modes that 
leave the horizon in the final N e-folds. This is well appre- 
ciated for scalar fields during inflation [32] ■ The modes 
that leave the horizon add up in a "random walk" man- 
ner to form a classical background that is experienced 
as homogeneous by modes of smaller size. A homoge- 
neous classical vector is a field that points in a given 
direction (that, in a given realization of the model, is 
determined by the random addition of the super-horizon 
modes) that breaks isotropy. In the magnetogenesis ap- 
plication, the effects of this energy (i?^) /2 on the back- 
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ground evolution have been well appreciated (see, among 
others, [26l HSl |44] ) . We point out that this energy is as- 
sociated with a classical vector field, and in this work we 
show that this vector imprints a strong anisotropy to the 
power spectrum and bispectrum of C in all the applica- 
tions of ([2|. 

Specifically, we show that the natural value of asso- 
ciated to these modes is ~ 0.1 (respectively, ~ 0.01), if 
inflation lasted about 50 e-folds (respectively, about 5 e- 
folds) more than the final ~ 60 e-folds necessary to gener- 
ate the CMB modes. Generic models of slow roll inflation 
are characterized by a much longer duration of inflation, 
and, therefore, embedding the mechanism ([2| in one of 
these models generically results in too anisotropic pertur- 
bations. For a tuned duration of inflation the mechanism 
becomes extremely predictive since, there is essentially no 
free parameter in ([2|. The only relevant quantity is the 
magnitude of the classical vector field present when the 
CMB modes left the horizon, and that can be "traded" 
for (7*. Therefore, any given value of should be asso- 
ciated with firm predictions for other observables. 

There are two such predictions that immediately come 
to mind: the first is a TB and an EB mixing in the CMB 
data |45| , resulting from the coupling between scalar and 
tensor modes that, due to the anisotropy, takes place al- 
ready in the linearized theory (TTJ H5H48j . The second 
is a directionality dependence in the bispectrum (and, 
in principle, in the higher point correlation functions), 
with a clear correlation with the one in the power spec- 
trum. The bispectrum resulting from ([2]) is computed for 
the first time in the present work. |^ We show that the 
isotropic power spectrum and bispectrum of |39j are in 
fact the theoretical expectation (i.e., the theoretical av- 
erage over several realizations) for the anisotropic signals 
that we obtain here, and which are the real quantities 
that are produced by any single realization of Quite 
interestingly, an observable g^, produced from ([2j) is as- 
sociated to an observable bispectrum which is enhanced 
like the local one in the squeezed limit, and which has a 
characteristic shape and anisotropy (immediately corre- 
lated with the one in the power spectrum). 

The paper is organized as follows. In Section [ll] we 
study the spectrum of the vector field perturbations ob- 
tained from ([2| , with a particular attention for the func- 
tions / that result in scale invariant vector modes. In 
this section we also further discuss the role of the large- 
wavelength modes in determining the classical back- 
ground anisotropy that affects modes of CMB wave- 
lengths. In Section |III| we study how these modes 5A 
are coupled to the modes of C, through ([2]). The power 
spectrum and bispectrum of C, are computed respectively 
in Sections |IV| and [V] There, we show explicitly how the 
sum of the long- wavelength modes adds up with the solu- 
tion of the classical equations of motion to determine the 



^ For previous works on anisotropic non-gaussianity, see |284432l 

IMll]. 



physical value of ^'classical (or ^classical) observed by the 
CMB modes. The resulting phenomenology is reviewed 
in Section [VTl In Section [Vll| wc discuss our results, that 
generally apply to all the realizations of ([2| that give 
a nearly scale invariant vector field, in the context of 
anisotropic inflation [19] , of the waterfall mechanism [28] , 
and of magnetogenesis [251 HO] ■ A concluding discussion 
is given in Section [VlII| 



II. A SCALE INVARIANT VECTOR FIELD 

Let us consider a locally U(l) invariant vector field 
with lagrangian ([2]). Ref. [35] identified this field with 
the electromagnetic one, assuming that / sets to a con- 
stant after inflation (in this case, we can simply normalize 
^ond = !)■ The function I enters in the definition of the 
electric and magnetic components 



E, 



(I) 



A' , B, 



(I) 



£ijkdjAk 



(3) 



(we denote by (...) the vacuum expectation value of a 
field, or of a function), where prime denotes derivative 
with respect to conformal time r, and a is the scale fac- 
tor of the universe, ds^ = (r) (— dr^ -I- dx^). With the 
notation the physical energy density in the vector 

field assumes the conventional expression p — ' ' ^' — - 
at all times. In this work, apart from where we explic- 
itly refer to the magnetogenesis application [25] , we do 
not necessarily identify the vector field with our pho- 
ton, but we keep the "electromagnetic" notation ^ for 
convenience. The classical equations of motion obtained 
from ([2]) are solved by a homogeneous "electric" field 



oc 



^j2-py. While the standard case, / = const, cor- 
responds to pe oc a^^, a constant "electric" energy is 
obtained if (/) cx a^^. 

A desired time evolution for (/) can be obtained for 
several functions I {<p), provided they are suitably ar- 
ranged with the infiaton potential [5(7. Indeed, 



a (X exp 



dip 



where we have introduced the slow roll parameter 



p 



(4) 



(5) 



(prime on a function here denotes derivative with respect 
to its argument) and where a monotonic slow roll inflaton 
evolution with < is assumed. Therefore, a desired 
behavior (/) = / (a) can be obtained by choosing the 
functional form of / to coincide with that function / of 
the right hand side of Q: 



I = Iq exp 



n dip 



iv>)Mp 



(I) 



(6) 
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where Iq can be chosen so that 7=1 after inflation. As 
a concrete example, the choice 



/ = e 



2M^ 



(7) 



results in (/) cx a^^^. We thus see that a constant i?'"-' 
is achieved for c = 1 . Many other choices of V and / are 
clearly possible, provided that their functional forms are 
related to each other to produce through Q the desired 
time dependence for (/) [iP] . 



A. Production of vector fluctuations from the I F 

term 

It is well known that, for / = 1, the vector field is 
conformally coupled to a FRW background, and so its 
fluctuations are not excited by the expansion of the uni- 
verse. On the contrary, as we now review, the choice 
(/) oc that allows for a constant i^'-"-' solution, also 
excites the vector fluctuations to produce a classical and 
scale invariant spectrum of "electric" fluctuations at large 
scales. 

We quantize the vector fleld in the Coulomb gauge 

= 



V 



E 



ax (fc) Vx (fc) + 4 (-fc) (fc) 



(!) 



(8) 



where ex are circular polarization vectors satisfying the 
relations k ■ e± (jej = 0, k x e± (jtj — ^ike± (jtj , 

e± (^^^ — ^± (j^^ J ^^'^ normalized according to 
ex (^kj ■ ex' (jtj = Sxx' ■ The annihilation / creation op- 



= Sxx' S'^''^ k-k 



erators satisfy ax yt 

The mode functions satisfy the evolution equation 



{I) 



Vx = 



(9) 



where prime denotes derivative with respect to conformal 
time T. For (/) cx cx (we disregard slow roll 
corrections, so that a — —jj:p), the properly normalized 
vector modes are 



Vx 



1 -|- ikr 



-ikr 



(10) 



We Fourier transform the "electric" and "magnetic" 
fields ^ 



E = 



3/2 



B 



d^k 



{2n) 



3/2 



(11) 



Inserting the solutions ( 10 1 in ([s]) we see that the "elec- 
tric" and "magnetic" fields become classical (commuting) 
fields at super-horizon scales 



5E (fc) = ^fe (fc) [ax {k) + a 

X 

5B{Pj =Y.Bk \ex (^) 

A 

3i72 



ax [kj +a\(-k 



£k ^ — 7- , Bk ^ — r- , -kr < 1 (12) 

Namely, the "electric" field fluctuations are nearly con- 
stant outside the horizon, while the "magnetic" field fiuc- 
tuations rapidly decrease. As we mentioned, the electric 
field fiuctuations are scale invariant (slow roll corrections 
will slightly tilt their spectrum; however, we disregard 
slow roll corrections in this work whenever compared 
with a non-vanishing expression at 0— th order in slow 
rofl). 

Finally, we note that this mechanism enjoys a du- 
ality symmetry (/) ■(-> Under this exchange, the 
"electric" and "magnetic" modes interchange their role. 



SE 



6B 



The original mechanism [5S] aims to 



produce fiuctuations with scale invariant magnetic en- 
ergy, and therefore has (/) cx during inflation. This 
corresponds to choosing c = — 1 in the example ([7|. 
For definiteness, our explicit computations are done for 



(I) 



However, our results can be readily extended 



to the context of [25] by exploiting this duality. 



B. Classical anisotropy when the CMB modes 
leave the horizon 

We denote by the "electric" field obtained from 
solving the classical equations of motion of a given model. 
For instance, as we discuss in Subsection |VII A the clas- 
sical equations of motion of the anisotropic infiationary 
model [in], characterized by c ~ 1 in ([^, admit an at- 
tractor solution with a nearly constant E'-^^ . 

We denote by i?ciassicai the classical and homogeneous 
electric field measured by a local observer at some time t 
during infiation. We can assume that at the initial time 
of infiation Ti„ no classical fiuctuations are present, so 
that i^ciassicai = E^^"^ at Tin. Howcvcr, this identification 
is no longer exactly true at any later time. Indeed, at the 
time r > during infiation 



E, 



lassical 



+ E 



IR 



(13) 



where the second quantity (IR = "infra-red" ) denotes the 
sum of all the modes SEk that left the horizon between 
the times Tin and t. These modes have become classical 
and are homogeneous from the point of view of a local 
observer present at t. 

The same considerations apply to the "magnetic" com- 
ponent of the vector field. For any single realization of 
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the mechanism (|2| with (/) cx a ^, the quantities E^^ 
and are drawn by, respectively, a gaussian (to very 
good approximation) statistics, with vanishing mean and 
with variance 



2 



{SE {xf 



1 



dk\£kf 



9H* 



IR 



dk 



27r2 ./jR k 



9H^ 
2^ 



TV 



2 _ ■3ir_ 



(14) 



where the IR modes are those characterized by momen- 
tum k in the interval < fc < — , and where N is the 

— Tin —T 

number of e-folds of inflation from Tin and r. 

These quantities are the natural expectation for the 
energy that gets progressively stored in the "electromag- 
netic" field during inflation 



PSEM 



{SE (f ) 



47^ 



N 



PB < PSE 



(15) 



This has been well appreciated in the magnetogenesis ap- 
plications of this mechanism (we remark that the role of 
the "electric" and "magnetic" energy is interchanged for 
(J) ^ -j^). For example, ref. computed the energy 
density accumulated in these super-horizon modes for all 
possible values of n in the (/) oc a" dependence. For 
\n\ > 2, this energy density grows as a'^d^/^l-i) (while 
n = ±2 results in p oc Ina = iV, as we have seen), until 
the backreaction of this energy is no longer negligible. 
Ref. HHj studied the regime of strong backreaction that 
takes place for n > 2. Ref. [Hj studied magnetogen- 
esis for a more general time dependence of /, imposing 
as one of the conditions that the energy density in the 
classical super-horizon modes remains less than that of 
the inflaton. 

We obtain a first limit on the total duration of infia- 
tion A^tot by imposing that E^^ has a negligible effect 
on the background evolution. The strongest backreac- 
tion constraint does not actually come from psE <SC 
but rather from the evolution equation of the inflaton 
(loosely speaking, it is "easier" for the vector field to af- 
fect the motion of the inflaton, that is slowly rolling on a 
flat potential, than the expansion rate). The correspond- 
ing condition is psE ^ 2 e , and it is satisfled for 
TVtot « 0(10^) [39 . 

The variances of the three components of E^^ are equal 
to each other, |^ and equal to 1/3 of the value given in 
(14). This, however, does not mean that in a given real- 



in all three directions. Indeed, the difference between dif- 
ferent directions is drawn from the statistics 

{SEl - SE'^) = 

y((<5i?2 - SE^^Y) ^ 2 (SEl) = ^ psE.N (16) 

Therefore, psE is also the typical amount of the classical 
anisotropy provided by the IR modes. This is not sur- 
prising since, in any given realization, E^^ is a classical 
vector that points in some given direction. 

Only for > \E^^\, or, equivalently, for p^io) » 

PSe , the classical electric field measured by the observer 
is (deterministically) given by the solution of the classical 
equations of motion. If this is not the case, one should 
conclude that the solution of the classical equations of 
motion is unstable under quantum corrections, and one 
should expect large corrections to the predictions made 
if only is considered In our computations below, 

both the contributions to -^classical will be accounted for 
(as we will see, this amounts in considering loop contri- 
butions to the power spectrum and the anisotropic spec- 
trum of the cosmological perturbations). 



III. ANISOTROPIC SOURCE OF THE 
COSMOLOGICAL PERTURBATIONS 

As we have discussed in length in the previous Sec- 
tion, for the mechanism we are studing the cosmological 
perturbations that we observed experienced a classical 
homogeneous vector field (13) when they left the hori- 
zon. This breaks the background isotropy, and, strictly 
speaking, the local patch where these modes live has a 
Bianchi-I geometry with residual 2d isotropy 



-dt^ + (t) dx^ + (t) [dj/ + dz^ 



(17) 



where for the definiteness the a;— axis has been oriented 

along i?classical- 

To characterize the anisotropy, we define 



If 



3 {Hy — H'X 



(18) 



This nearly corresponds to the mechanism of anisotropic 
inflation of [TS], where the model ([t]) has been employed, 
and where it is shown that the classical equations of the 
model admit an anisotropic solution with ~ (c — 1) e. 



ization the super-horizon modes add up to equal amounts 



In [19 , and in the successive works j34H57j that study the 
linearized perturbations of this model, it is assumed that 
-^classical — E'^^\ while, as we have discussed, these two 



* Actually, one can obtain a small difference proportional to the 
small difference AH in the expansion rates of the different direc- 
tions. As we well see, AH/H needs to be <J 10~*, and therefore 
this difference is completely negligible for the present discussion. 



^ Before eq. | |l3| , we set -Eclassical = E^^' at the start of inflation. 
We note that if this is not the case the departure of the classical 
electric field from the solution of the classical equation of motion 
will in general be even greater. 
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quantities are in general different. Tlierefore, this model 
leads to anisotropic inflation even for c = 1. Given this 
strong correspondence, our study has several relations 
with 519], and for example we will show that our results 
for the power spectrum coincides with that of 35 37J, 
once the value of E'^^^ used there is replaced by the full 
^-classical (for the bispcctrum instead, our result is com- 
pletely new, since the perturbations of |19j have been so 
far studied only at the linearized level). 

In the studies [51H57| , the perturbations are separated 
according to how they transform under an S0(2) rota- 
tion in the symmetry plane orthogonal to the vector vev, 
as originally done in [TTJ UT] . This simplifies the prob- 
lem, as modes that transform differently are not coupled 
to each other at the linearized level. Still, the non-FRW 
background results in a very involved computation, once 
the perturbations of all the fields (metric included) are 
taken into account. In particular, the anisotropy does 
not increase the number of physical modes, but results 
in couplings between these modes that would be absent in 
the FRW case [HIIIT]. The main result of [55^57]. which 
is somewhat surprising a-priori, is that an anisotropic pa- 
rameter 5* = O(O.l) is obtained for a ^ = O (lO^®) 
background anisotropy. Motivated by this result, one can 
instead use a different approach, and perform the stan- 
dard quantization of the cosmological perturbations on a 
FRW background, ignoring at zeroth order the couplings 
between the different modes. Such couplings can be taken 
into account as perturbative mass insertions in the in-in 
formalism. This is effectively the procedure adopted in 
[35. ,37j when they solve analytically the linearized theory. 
Moreover ref. [37] showed that the dominant operator 
that determines g^, in this perturbative evaluation is the 
dip—SAf^ coupling obtained from expanding the vector ki- 
netic term, and with no contribution from the metric per- 
turbations. This analytic approximated result is in excel- 
lent agreement |37| with the one obtained from an exact 
numerical evolution of (in which the full quadratic ac- 
tion of all the 2d scalar modes is retained) . Moreover, the 
analytical and numerical results of [371 agree with those 
of [3S] and respectively. 

We use the same computational scheme for the bis- 
pectrum computation. Specifically, we disregard metric 
perturbations, and we use in the in-in formalism the 0th 
order eigenmodes obtained from the approximate FRW 
quantization. The bispectrum is produced by interac- 
tions, and we know that for slow roll FRW inflation the 
interactions of the metric perturbations produce an un- 
observable bispectrum. We will see instead that the in- 
teraction between the vector and the scalar field, which 
is encoded in the vector kinetic term, results in a larger, 
and potentially observable, signal. Ref. [39 proved this 
explicitly for the case of -^classical = 0, by solving the 
second order equation for the curvature perturbation in 
spatially flat gauge ( — —^6ip (we recall that prime de- 
notes derivative with respect to conformal time r, while 
= ^). It was shown in [35] that the contribution 



from the direct vector-scalar interaction is slow-roll en- 
hanced with respect to that coming from the interactions 
of the metric. There is no reason to expect that the rela- 
tive strength of the effects should change for a O (10~^) 
background anisotropy. 

Therefore we retain the FRW quantization of the vec- 
tor field performed in the previous Section. The remain- 
ing perturbations (of the scalar field and of the metric) 
are also quantized as in the standard FRW case. We 
then expand the vector kinetic term into SA and Sip. We 
use the resulting interactions to determine the dominant 
anisotropic contribution to the power spectrum of ^ and 
the dominant contribution to the bispectrum. 

Therefore, all the dominant effects arise from expand- 
ing the only interaction term between the infiaton and 
the vector field, 

(19) 

In spatially flat gauge, 5(p = C — ^/2eMpC. We 
note that, in principle, C has additional contributions 
proportional to the perturbations of the vector field. Ref. 
Ii9 showed that these contributions are completely sub- 
dominant in the case of ii'ciassicai = 0. We believe that it 
is very natural to assume that this continues to be the 
case also in the current context. Indeed, as we shall see, 

^ needs to be < O (lO~®) during infiation or otherwise 

Pip \ / 

the power spectrum is too anisotropic. This ratio further 
decreases between the end of inflation and the inflaton 
decay, when the inflation field performs coherent oscilla- 
tions, so that Pip (X a~^, while the vector kinetic term 
becomes standard, and ps oc a^^. If this is the case. 




where T^ch and a^eh are, respectively, the temperature 
of the infiation decay products and the scale factor at 
the inflaton decay. We see that it is therefore natural 
to disregard pE and Sps at reheating. In this way, the 
only relevant contribution of SA to the final curvature 
perturbation is the modification of Sip induced by the 
vector-inflaton coupling (which is precisely the effect that 
we are computing). We note that this assumption is also 
made in |35H37| when they give the power spectrum of C 
in the model [19]. 

Using the expression ([6|, we have, for the first two 



We note that metric perturbations are also disregarded in the 
computations ^28 33 of the anisotropic bispectrum through the 
waterfall mechanism. 



(1) (1) 



+ 



+ 



/C1 
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FIG. 1. Leading diagrams for (C^), with the vertices labelled as in (26 1 



terms in the expansion of 



-2n 



6^^ 



6M2 



2eM„ 



AP 



Sip = -2n{P)C 



2e 



2n^ (P) C 



(21) 



where 77 = ^ is a slow roll parameter, and where in 
the final approximation we retained only the dominant 
term in slow roll approximation. 

We inserted these expressions in (19), taking n — —2 
(which corresponds to a constant "electric" field). We 
expanded also the vector part as in ([s]), and we com- 
puted the contributions to the power spectrum and the 
bispectrum combining the resulting vertices. We veri- 
fied that the leading diagrams do not contain interactions 
coming from the second order term (if this was not 
the case, we should worry about the convergence of the 
/ [(</?) + Sip] expansion). Therefore we have the two dom- 
inant cx SipA'^'^'f SA and oc SpSA^ interactions!^ 



£int 3 a"* 



AEi'^^E^C + '^SE-SEC^ 



(22) 



We note that the vector field enters in the quadratic 
term £int.i only in the combination cx 2{Ex) SE^ C 
— B'^ (X F^^Ff^". We also note that the cubic term 

'Cint,2 should also have the term —2a^5B ■ SB C; this term 
however gives a subdominant contribution to C with re- 
spect to contribution from the "electric" components, 
since the "electric" modes are much greater than the 
"magnetic" ones at super-horizon scales - see eq. ( 12 ) 



- which is where they provide the greatest contribution 
to C (as discussed in [35] and below). 

We stress that, although in this Section we have often 
mentioned the anisotropic infiationary model |19| . this 
interaction lagrangian applies to all models that verify 
the two assumptions spelled out at the beginning of the 
introduction. We actually see that the precise functional 



forms of / and V do not enter in (22), apart from the 
fact that the relation ([6| has been imposed, with n ~ —2. 
Identical results would be obtained for n = 2, exploiting 
the "electric" - "magnetic" duality of the mechanism. 



We recall that n = ±2 precisely correspond to enforcing 
the assumption 2. made at the beginning of the Intro- 
duction. For the choice ([7|), one has n = —2c. However, 
we stress that (22 1 is valid independently of this choice. 



The interaction lagrangian ( |22[ ) enters in the n— point 
correlation functions through the in-in formalism relation 



.C?'(r), i/i,,(ri) 



dr 



N 



Hint {tn) 

(23) 

where in our computational schemes the quantity C^"'' is 
the Fourier transform of the (unperturbed) FRW quan- 
tized field 



c 



(0) 



d^k ^(0) 



(2^) 



3/2 



;(o) 



^k 



- + a 



(0) * J 



At leading order in slow roll, we have 
(0) _ H{l + ikT) 



(24) 



(25) 



The unperturbed fields are also employed in the in- 
teraction hamiltonian Hmt (t) = — J d^xCmt {t, x). The 
two terms in ( 22 ) give, respectively, rise to the two terms 



i^int.l (t) 
Hint,2 (t) 



(0) 



d^kSE^ (r, Pj C^lii 



d'^kd^p 

{2TTf'^ 



5E 



(r,fc).5^(r,p)Cil.(r) 
(26) 



Once inserted into (23), the two interaction terms give 



rise to the leading contributions to the power spectrum 
and bispectrum described by the diagrams shown, re- 
spectively, in Figures ([T]) and These quantities are 
computed in the following Section. 



IV. ANISOTROPIC POWER SPECTRUM 



When we perform the expansion, linear terms in the perturba- 



tions are removed by the background equations of motion; we 



8 



The total power spectrum I fc j is related to the full 



two-point correlation function ( 23 ) by 

6^^^ (fci + h) 



{ki) (27) 



We denote the contributions of the first and of the last 
two diagrams in Figure [T] as, respectively, 



and, correspondingly, 



SP 



(28) 



(29) 



These quantities are computed in the following Sub- 
sections. 

A. Tree level contributions 

The first diagram in Figure [l] gives the (unperturbed) 
FRW power spectrum 



■p(o) ^ 



ij4 



87r2eM2 



(30) 



at super-horizon scales, where the slow roll expression 
(25 1 has been used. 

The second diagram in Figure [T] gives the anisotropic 
contribution 

*(4,4(t»|i = - /' *i ' tin 



, "Hint,! (^2) 



(31) 



We are interested in the power spectrum at super horizon 
scales, fc|r| <C 1. The time integral ( [3l] ) is dominated by 
the times for which the modes in Hi^t (^i) are also out- 
side the horizon (mathematically, the contribution in the 
sub-horizon phase is suppressed by the oscillatory phases 
in the mode functions). This condition will be relevant 
for setting Tmin (see below). As remarked before eq. ( 12 1, 



in this regime the vector field is classical, and does not 
contribute to the commutators. The only nontrivial ele- 
ments in the commutators are therefore 



Cf (r) , (r') 



:(0) 



(Cf (r)Cr(r')-c.c.) 



fc' 



-iH 



6eM2 



(32) 



where the last result is true in the super-horizon regime. 
We insert (|26l) into (31) and perform the commutators 



between the fields. The two resulting (5— functions 
are employed to perform the two integrals over momenta, 
leading to 



(0)2 



^ i—1 



n 



~Z4 



5E, n, kA5E., T2, fc2 ) (33) 



Requiring that the vector field in this expression is in the 
super-horizon regime limits each time integral to > 
This sets the value of Tmin in the two integrals. 



ki ■ 



Using the expressions (12), and the identity 
- (^) 



A 



(^) 



ki kn 



(34) 



we obtain 



^(4,4(-))ii 



,(0)2 5(3) (^k^ + fc2^ 



2E.. 



kl 



sm 





"^3 _ ^'3 








Changing variable y' = and recalling that — fcir <C 
1, the time integral in the second line becomes 



I. 



dy' 
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In- 



y 



— kl T 



(36) 



At the end of inflation this quantity becomes N-^^ , namely 
the number of e-folds before the end of inflation at which 
the modes with wavenumber k\ left the horizon. Using 
this result. 



(37) 



where Tend denotes the end of inflation (we assume that 
/ rapidly approaches 1 at the end of inflation, so that the 
vector field is rapidly diluted away by the expansion, and 
the power spectrum of C freezes out). 



For the anisotropic infiationary model of [TH] , eq. ( 37 ) 
coincides with the analytic result of [33 [13 ' which, as we 
remarked, is in excellent agreement with the full numeri- 
cal computation of [ 36. i37j . This confirms the validity of 
all the approximations that we have performed. 



B. Including the Loop contribution 



commit a mistake by disregarding the effect of the gauge field 
vev in the background evolution, but, as we remarked, we work 
in a regime where this effect is negligible. 



The expression for the loop diagram in Figure [Tj is 
analogous to eq. (31), with "Hint, 2 replacing ?^int,i- We 



perform the commutators, again keeping in mind that 



9 



the time integrals are dominated by fields in the super- 
horizon regime. We obtain 



1 



5Ei (ti, p) SEi (ti, ki ~ pj SEj (t2, q) 6Ej (t2, h - g) 



(38) 



Evaluating the full correlator (disregarding the discon- 
nected contraction) gives 



theory 



2e2Af4 



{2Trf p^\ki-p\3 



n 



dTi 



(39) 



where also in this case we recall that the main con- 
tribution comes from the times when all the fields are 
in the superhorizon regime (for the same reasons men- 
tioned after eq. (31)). This corresponds to T^-m — 

Max ^ 1 



p' |fei-p|. 

keeping only the logarithmically enhanced term, leads to 



Performing the time integrals, and 
mically enhanced t( 



theory 
2 , 



2e2M4 



W3 1 + cos2 6'.t: ^ / 



1 



1 



-rp -r|fci-p1 



(40) 



The final momentum integral naively diverges logarith- 
mically at the two poles p — ?► 0, k\. We need however to 
impose that the the fields 8E in ( 38 1 were inside the 



horizon at the start of inflation (or they would not be 
excited by the mechanism described in the previous Sec- 
tion). Therefore p, |fci— p| > where Tin denotes the 
conformal time at the start of inffation. With this cut-off 

bP2 (rend, fc) Ithcory ^ 192 V'^''^'^ ^1 {Nt^t " Nk) (41) 

This result accounts only for the contribution of the two 
poles and therefore is accurate only as long as the log- 
arithmic enhancement encoded in the last term is 3> 1. 
For A^tot — this last factor should be replaced by a 
0(1) factor. The result (41 1 was first derived in P9] 



using the Green function method. 

We note that the modes contributing to the logarith- 
mic enhancement are precisely the modes that left the 
horizon during the first TVtot — ^cmb e-folds of inflation, 
and that contribute to the classical field E^^. As we 
discussed after eq. (13), this quantity adds up with the 



solution of the classical equations of motion to give 



-^classical, which is the homogeneous and classical "elec- 
tric" field observed by the CMB modes. The quantity 
( 41 ) is the theoretical expectation value associated with 



the loop diagram; it is the result that one would obtain if 
one could average over several realizations of the mecha- 
nism (see also |56j). However (assuming that this mech- 
anism describes our universe) the CMB modes exit the 
horizon after a single realization of the first iVtot — -^cmb 
e-folds of inflation. They therefore are not affected by 
the theoretical average (|i?^^| ^) (which is statistically 
isotropic), but by the value that E^^ happened to as- 
sume in that single realization. 

Therefore, if we want to compute the contribution of 
dP2 to the power spectrum that is observed in a sin- 
gle realization, we need to replace the quantum opera- 
tor of the IR modes entering in ( 38 1 with the classical 
Fourier transform of E^^. If we do so, the expression 



( 38 ) becomes formally identical to (|33| , with the quantity 
^;(o)2 ^ 1^(0) |2 replaced by \E^^f and with 6E^ (r^, k,^ 
(namely, the component of the fluctuations along the di- 
rection of E^'^'>) replaced by the component of SE ^r^, ki^ 
in the direction of |i?^^|2. As a consequence, the value of 
'5^2 (''end J kj for a single realization coincides with (37 1, 
with replaced by E^^ , 



/ \ 24 li?^^|2 

SP2 (rend, kj I localization ^ " ^1^'"^ siu^ 

Computing the theoretical expectation for this con- 
tribution amounts in replacing sin^ 6f^ E^^^ with 

I a|,jp ^. The resulting expression coincides with ( j4T| . 

We note that taking the classical IR value for one prop- 
agator is equivalent to taking the classical IR value for 
one of the two vector fields in the interaction hamilto- 
nian i?int,2 given in eq. (26 1. This interaction term 
then becomes identical to i?int,i apart from the fact that 
is replaced by E^^. The two expressions than be- 
come a unique vertex in terms of the vector i?ciassicai = 
^(0) This is the expected physical result, since the 

CMB modes "measure" the sum i^ciassicai and cannot dis- 
tinguish the two components. Even at the diagrammatic 
level, taking the IR limit in the third diagram of Figure 
[T] amounts in shrinking to zero one of the propagators, 
and one thus recovers the second diagram with a different 
external vector in the vertex (clearly, all this discussion 
applies also to the bispectrum computation). Therefore, 
from the last two diagrams of Figure [T] we obtain 



5P (Tend, fc) ~ ^ 



— ^^^^f'^^' NiV^°^ sin^ 9f^ E^^^^^^,,^ (43) 



Viip) 



ANISOTROPIC BISPECTRUM 
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FIG. 2. Leading diagrams for (C ), with the vertices labelled as in (26 1 



The bispectrum of C is defined as 



Be (r, ku k,, h) 6^'^ (^Y. = {k k k (^) 

(44) 

The diagrams shown in Figure [2] give the dominant 
contributions to the bispectrum. The computation is pre- 
sented in the foUowing Subsections. 

A. Tree level contributions 



We evaluate the contribution of the first diagram Bi 
as we did for SPi in the previous Subsection. We start 
from 



i^.ki^Jr))^^^Jdr,dr,dr,{[l\C^C|^C|^(r), 

^int,2 (ri) ] , Hint,l (rz) ] , Hint.l (T3) ] ) + . -(45) 



where the dots denote two additional terms obtained by 
permuting the position of "Hint, 2- Performing the com- 
mutators between the fields results into 



^fcl ^fc2 ^fc,^ 



n 



1 27e^mAB 

^ i—l 



SEr, (ti, fci SEx (t2, fc2 



d^p 

SE{T3,p)-dE(T3,h-p) +■■■) (46) 



where dots denote two additional terms obtained by per- 
muting k^ with the other two momenta. As for the di- 
agrams evaluated in the previous Subsection, the time 
integrals are dominated by the times for which the mode 
functions are classical and do not oscillate. For the term 



that is explicitly written in ( 46 1 this gives the three lower 
limits Ti > ^ 



ki ' 



7-2 > 



and T3 > Max 



j_ i_ 

'fci' fe2 



(the limit on is most easily seen after taking the ex- 
pectation value). Taking the expectation value and per- 



forming the time integral, we obtain 

' Nk.Nk, Min [Nk,, Nk. 



Bl I Tend , fc, 



75(0)2 



1,3 jL3 
"■1 '*'2 



where we have defined 



■C 



fci k2 E(o'> ~^ ^ permutations 
(47) 



Ck„k2,V = 1 - COs2%^_^ - COS^ Of^^ y 

+ COS y COS %^ _ y COS _ (48) 



B. Including the loop contribution 

For the second diagram in Figure [2] we obtain 

Cfci Cfc^ Cfc3 (T-cnd)^ 
' 3 



2 27e^H^M^ 



n 



,3/2 



dr, 



-SE {Ti,pi) ■ SE (Ti,ki -p, 



(49) 



where again we restrict the time integrals to when the 
modes are outside the horizon. The full correlator gives 



^fci ''£2 ^£3 V^cnd j Ithcory — 



(2^) 



9/2 



JL f -T^ Qki [p] Qij 

t—1 '■ 



P~ki 



p3 Ip*- fci|3 |p-f /C3 



3 



(50) 



where we have defined Qij [p\ 



Pi pj and 



where the time integrals are restricted to ri > 



Max 



, T2 > Max 



|p-fci|' Ip+fcal 



and 



T3 > Max , , ^ , 

The final momentum integral has three poles where 
it naively diverges logarithmically. The integral is regu- 
lated by the same argument used after ( 40 ) , which in the 
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present case enforces p, \p — ki\, |p + fcs | > 
forming the integrals gives 



Per- 



^2 Tend, 



i j I theory 



1152y27r5/2-p(0)3|^^^ 



xMin [Nk, {Ntot -NkJ, Nk, (A^tot - Nk,)] 



+2 permutations I 



1 + COS^ I 
1.31.3 



felfc2 



spectrum. Conversely, this indicates that the classical 
"electric" field must be extremely subdominant not to 
conflict with the phenomenological limits on the isotropy 
of the power spectrum. 

We recall that, in any given realization, i?ciassicai = 
^(0) _|_ j^iR^ where E^^^ is the solution of the classical 
equations of motion of the model, and is drawn by 



(51) 8- gaussian statistics of zero mean, and variance (14) 



As for (41), this result is accurate for iVtot N/^.; if 
this is not the case, the factor iVtot — should be re- 



placed by a O (1) factor. The result (51 ) was first derived 



in [39^ using the Green function method. 

As for SP2, the result is dominated by the regime in 
which one propagator has an IR mode. Also in this case, 
the value for a single realization is obtained by replacing 
the quantum operator of the IR field entering in ( 49 ) by 



a classical Fourier tra nsfo rm. We recover an expression 
formally identical to (46), with E^'^'> replaced by E^^ . 



Therefore, precisely as for the power spectrum, the total 
contribution of the two diagrams to the bispectrum is 
given by 



B I Tend , ki 



288\/2 7r5/2 I ^classical I 



■p(0)2 



^ Nk.Nk, Min [Nk„ Nk,] 

^ ^ h3 1,3 ^*:i,fe2,-E=i='=='=^»i 

'>'l '^2 



2 permut 

(52) 



where we recall that C is given in ( 48 ) 



VI. PHENOMENOLOGY 



The total power spectrum is given by 



1 



94 fT'^ 

classical 



v{^) 

corresponding to a negative 



94 I7'2 

_ classical 7X7-2 

94 

^ classical ]\t2 




parameter in ([T]) 



94 

classical 



(53) 



(54) 



where the approximation is due to the fact that \g^\ is 
phenomenologically constrained to be << 1. 

Inverting this relation for the CMB modes, and denot- 



ing by pE,, 



the energy of the classical "elec- 



tric" field present when the CMB modes left the horizon 



5.8- 10 



-9 



\9- 



*|CMB 



0.01 0.1 



60 



N, 



CMB 



(55) 



and therefore we see that even a very subdominant vec- 
tor field can produce an appreciably anisotropic power 



E 



IR 2 



N 



2^ 



N 



(56) 



The variance grows with the number of e-folds of infla- 
tion. 

It is instructive to evaluate the ratio 



Rn 



^IR 2\ 

/ N 



(57) 



PSE,N 

at the moment that the CMB fluctuations left the hori- 

37 



zon. Combining (55) and (56), and using the observed 
value [5 for the dominant term V^^^ ~ P^- ~ 2.5 • 10~^ 
we obtain 



R 



|g*|CMB 

0.1 



60 

A^CMB 



CMB 



(58) 



The numerator of (57) is the observed value, while the 



denominator the theoretical variance. Therefore, a value 
RNtot-NcMB — 1 should be naturally expected for this 
ratio. Indeed, R^t^^-N^f^^ ^ 1 indicates either that the 
sum of the IR modes is unnaturally small in that realiza- 
tion (i.e. 1^(0) |2, |i?iR|2 < /^IR 2\ )^ ^j^^^ 

\ ^ / Ntat-NcMB 

it unnaturally cancels against . Assuming a natural 
-RATtot-ATcMB ^ 1 reahzation, 



Is* I CMB Min 



0.1 



N, 



CMB 

60 



37 



CMB 



(59) 



In generic slow roll inflationary potentials , the total du- 
ration of inflation exceeds (by much) the minimal A^cmb 
amount. One naturally finds an order one anisotropy in 
these models. 



Another important conclusion can be drawn from ( 58 ) . 
As it is exponentially unlikely to have jS^^p much 
greater than its variance, a value RNtot-Ncm-B ^ 1 i™" 
plies that lE^^^ is much greater than the sum of the 
IR modes. Only in this case the solution of the clas- 
sical equations of motion provides an accurate value 
for the classical vector field that affects the observables 



* Strictly speaking, eq. \55\ , and those derived from it, are only 
valid for \gt\ < 1. A value of A^tot — Acmb that results in a 



> 1 in ([58| should be associated to a 
natural -RjVtot-JVcMB ^ ^ regime. 



0(1) in the 
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curvature perturbations (as it is typically assumed in 



many works). A value R 



Ntat-Nci. 



1 indicates in- 



stead that the contribution of the IR modes is no longer 
negligible, and that the solution of the classical equa- 
tions of motion is unstable under quantum fluctuations 
(more appropriately, under the sum of the quantum 
fluctuations that have become classical). We see from 
(581 that RNtat-NcMB ^ 1 be obtained only for 
iVtot - NcMB < 37. 

Therefore, this mechanism can result in a |(7*|cmb — 
0.1 anisotropy for either a tuned short duration of zin- 
flation or an unnaturally small i?ciassicai- Assuming that 
this is the case, we can obtain a firm prediction from 
the anisotropy, the shape, and the magnitude of the 
non-gaussianity in the model. The bispectrum from this 
mechanism is given in eq. ( 52 ) . To quantify whether this 



bispectrum is of observable magnitude, we notice that it 
is enhanced in the squeezed limit precisely as the local 



template. We therefore insert (52) into the relation 



Be 



10 



{27: /nlPc 



(60) 



which, for the local bispectrum template, results in a con- 
stant /^L^' with the correct normalization (the non con- 
ventional factor of "27r-dependence" is due to our choice 
(24 1 for the Fourier transform). In the squeezed limit. 



we obtain 



/NL 



480 P(0)2 pE^, 



ki,k2,V 



fci < fc2 ~ 

, , . , , . . (61) 

where C is given in im. Recalling that V dominates 



the power spectrum, using (55 1, and simply denoting 
= Nqmb (assuming that the momenta are 



not too hierarchical; this is certainly the case in the CMB 
analysis), we arrive to the final estimate 



/] 



NL 



26 



ff*|CMB A^CMB '■^ki,k2,V 



Ct 



0.1 



60 



4/9 



fci < fc2 — fcs 



(62) 

where the value 4/9 has been inserted in the last factor 
so that this factor averages to one if we average over all 
directions of ki and k2 — k^. This result holds indepen- 
dently of the type of inflationary potential considered. 
In the estimate (62) we have averaged over all directions. 



with the assumption that this is what is done when ex- 
tracting the value of /nl from the full-sky data. The esti- 
mate ([62]) indicates that a |(;*| in the 0.01 — 0.1 range from 
this mechanism can likely be associated to a detectable 
bispectrum. A detection of a non- vanishing /nl of the 
local shape will motivate a more detailed analysis, where 
the full shape and the anisotropy of (52) are retained. 



VII. APPLICATIONS 



fields are employed to either break statistical isotropy or 
to generate primordial magnetic fields. 



A. Anisotropic inflation 

Ref. fT^ studied anisotropic infiation from the model 
^. It was shown that the classical equation of the sys- 
tem admit an attractor solution characterized by the ge- 
ometry (17 1 and by 



This values corresponds to [19 



AH 

If 



(63) 



(64) 



The power spectrum in this model was computed in 
[35H37] that however disregarded the contribution of E^^. 
Our computation reproduces their result (which confirms 
the validity of our computational scheme and the ac- 
curacy of our approximations) in this limit. Moreover, 
we have computed for the first time the bispectrum of 
this model, using the same computational scheme used 
to compute the power spectrum. 

The contribution of .B*^"-* to the anisotropy of the power 
spectrum forces c — 1 < 10~^ [551 - 157] . The implicit as- 
sumption in this result is that choosing a c — 1 — > in the 
lagrangian would result in a — >■ parameter. However, 
taking also E^^ into consideration, leads to the general 
conclusions on the natural amount of and on the in- 
stability of ( 63 ) that we have presented in the previous 
section. 



B. Waterfall mechanism 

Ref. [28] embedded the mechanism ([2| in hybrid infla- 
tion, through the lagrangian 

^ = ^ 2 ^ ^ 2 ^ 

^ 2 2 2,^ 2 2 , ^ u2 An A 2 

^9 X +^rn ip + -h M A^x 



(65) 



The additional coupling of the vector to the wa- 
terfall field X provides a contribution to C through the 
mechanism of modulated perturbations [57j . Specifically, 
the value of the vector field contributes to determine the 
end of inflation, that happens when 



\ 2 I 2 2 

— —Aw + g if 



h^A^A^' = 



(66) 



so that perturbations 5A^ get converted into C, at the end 
of inflation. This results in the additional contributions 
[28] 



In the following subsections we discuss how our find- 
ings apply to three well studied models in which vector 



A/] 



NL 



A5, -A^n (67) 
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where ip^ and T]e are, respectively, the value of the inflaton 
and of the slow roll parameter r\ at the end of inflation, 
and where we have normalized the coupling in ([2]) to 
(/) = 1 at the end of inflationj^ 

These additional contributions to the anisotropy ob- 
viously vanish in the limit in which the vector field is 
not coupled to the waterfall field (ft, 0). On the 
contrary, (59) and (62) are general results unavoidably 
present when the function / in ([65| is arranged to pro- 



vide a constant vector field (which is also necessary to 
have non- negligible values in (67)). These unavoidable 



contributions have not been included in the many works 
that studied [2B]. Given that ([59|) nearly saturates the 



phenomenological limit, the contribution (67) to is at 
most comparable to ( |59l). O n the contrary, the additional 
coupling ft present in ( 65 ) can allow for a greater non- 
gaussianity than (62). Indeed, in the minimal model ([2 ) , 
the parameter /nl is uniquely related to through ( 62 ) , 



and this strict relation increases the predictivity of the 
mechanism. For the model ( |65[ ), non-gaussianity can be 
enhanced for h ^ g. The shape of the contribution to 
the bispectrum resulting from the waterfall coupling |28| 
is in general different from the shape of (52); the two 
shapes coincide for h ^ g. 



C. Magnetogenesis 

Using the "electric" O "magnetic" duality that we 
have mentioned [57], the results for the perturbations 
that we have obtained are also valid in the case of 
(J) oc a^, which results in a scale invariant "magnetic" 
field [53]. In fact, the two results (41) and (51) have 
been first obtained in 



[39 for (/) (xcP. Although the 
mechanism ([2| was first introduced with this motivation, 
the magnetogenesis application suffers a strong coupling 
problem [53] . Indeed, with ([2| the fine structure constant 
scales as a oc (/)~^. The choice (/) oc corresponds to 
a fast decreasing a oc during inflation. Assuming 
that a reaches the present value at the end of inflation 
implies a ^ 1 during inflation, and a quantum field the- 
ory of electromagnetism that is (at the very least) out 
of computational control [26]. Alternatively, for a < 1 



at the start of inflation, one obtains a 



at the 



end of inflation, and requiring that a grows back to its 
present value before big-bang nucleosynthesis does not 
appear feasible [39]. 



Clearly, the strong coupling problem is a problem of 
the magnetogenesis application and not of the mechanism 
([2|, as one can always assume that a generic U(l) held is 
sufficiently weakly coupled to other fields. For a growing 
(/) the coupling decreases during inflation, and one can 
assume that the vector held has couplings < 1 at the 
beginning of inflation. For a decreasing (/) the coupling 
grows during inflation, and one can require a coupling 
<i 1 when inflation ends. In principle one may hope that 
a field with such properties is produced during inflation 
by this mechanism, and it is then (partially) converted to 
the magnetic one through some coupling. Some attempts 
so far in this direction have been unsuccessful |39j , but a 
general study remains to be done. If this, or some other 
idea "^T, can circumvent the strong coupling problem, 
our findings provide a signature that may be correlated 
with the magnetic field. 



VIII. CONCLUSIONS 

There are two questions associated to the attempt of 
reconciling a non vanishing anisotropy g^. with a model 
of inflation. The flrst one is "what are the other signa- 
tures that would accompany the measured value of g^, 
?"; the second one is "what is the natural value of 
in the model ?". A flrst difficulty that is encountered 
in answering these questions is that anisotropic hairs are 
typically erased by infiation [T5J [T3] , so that one has to 
design specific models that preserve the anisotropy. Vec- 
tor fields appear as the simplest possibility, as their vev 
breaks isotropy. Several mechanisms have been designed 
to prevent the rapidly erosion of the vector vev by the 
inflationary expansion. As shown in [18] . however, many 
of them (for instance, the use of a potential as in [14j . 
of a lagrange multiplier as in [S] , or of a coupling to the 
curvature as in [T5I and as in analogous models of vec- 
tor curvaton) have ghosts. The use of ([2| for preserving 
the vector vev [TU [28] has resulted in a rather excep- 
tional mechanism where the above two questions can be 
formulated. 

With the mechanism ([2|, the answer to the flrst ques- 
tion appears to be very positive. The anisotropy in 
the power spectrum is correlated with a characteristic 
and very likely detectable bispectrum. The bispectrum, 
whose full shape is given in ( |52| , is enhanced as the local 
one in the squeezed limit, where it has an effective local 
parameter (see ( 62 ) for more details) 



^ Ref. |32) studied the additional contributions to ( from tlic wa- 
terfall coupling in the case in which the vev of the vector field 
varies with time. For a discussion of the effects from a time 
varying inflaton coupling in modulated reheating see 1581 . 
^'^ A scale invariant magnetic field is also obtained at the start of 
inflation for (/) oc a~^. However, too much energy gets stored 
in the electric field in this case |26| . and the system enters in a 
strong backreaction regime in which a too small magnetic field 
is produced |43| . 



/l 



NL 



26 



|g*|CMB 

0.1 



(68) 



A larger bispectrum can be obtained if the vector held 
has additional interactions with the inflationary sector, 
as in the waterfall mechanism of [Ml- However these 



additional interactions are not needed, and the result ( 68 ) 



is the direct contribution that comes from the interaction 
([2]). Therefore, this is the result in the most minimal 
and predictive implementation of We see that an 
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anisotropy in the power spectrum at the 1% — 10% level, 
could be associated with a detectable bispectrum. 

For the second question, one needs to take into ac- 
count that, for a mechanism that results in a scale in- 
variant field (as in the implementations [IHIIIEDj a clas- 
sical background value of this field is unavoidably gener- 
ated by the sum of the modes that have left the horizon 
during inflation [42j . The crucial point for our discus- 
sion is that a background vector field breaks isotropy. 
It is unnatural to require that the anisotropy experi- 
enced by the CMB modes when they leave the hori- 
zon is smaller than the theoretical expectation value 

- (TVtot - Ncmb), where H is the Hubble rate 

during inflation, and Ntot—NcuB is the number of e-folds 
of inflation that took place before the CMB modes are 
generated, and during which the classical background of 
long-wavelength modes experienced by the CMB modes 
has accumulated. Taking this into account, we find that 
the natural value for the anisotropy is (see ( 59 ) for more 
details) 



Iff* I CMB ^ Min 



0.1 



CMB 



37 



(69) 



which can easily overcome the phenomenological bounds 
(particularly with the improvement expected from 
Planck [Slin])- We note that this saturation takes place 
for a very subdominant vector field, with very negligible 
backreaction on the background infiaton evolution. How- 
ever, as already obtained in |35H37| , even a very subdom- 
inant vector field can result in a substantial anisotropy 
parameter. 

Remarkably, these considerations apply also to the 



magnetogenesis applications of ^ . The role of the back- 
ground magnetic field has been so far disregarded in the 
study of the cosmological perturbations obtained in that 
context. However, the result (69 1 applies also to this 



case. We stress that the magnetogenesis application suf- 
fers from a serious strong coupling problem |26j . If this 
problem can be solved, a magnetic field through this 
mechanism would be naturally correlated with an ob- 
servable anisotropy of the perturbations. 

Analogously, an anisotropic signal should also be ex- 
pected in models where a triad of vectors is arranged 
to produce isotropic expansion or isotropic perturbations 
|59| . since the long- wavelength background values of the 
different vectors fields will in general be different. 

Although our study has been limited to the mechanism 
([2]), one may expect that an anisotropic background field 
and a large anisotropy of the perturbations may be a 
general outcome of all models that sustain higher than 
spin fields during infiation. For instance, it would be 
interesting to study the natural level of anisotropy to be 
expected in the mechanism of [22] (for which statistical 
isotropy has been obtained for very specific choices of 
the kinetic and the mass function [51]) or in the case of 
inflation or magnetic flelds from p-forms |60| . 
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